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Using the -light scattering matrix elements measured by a polar 
nepheiometer, a procedure for estimating the characteristics of 
ahsospherie particulates ms developed. A theoretical library data 
set of scattering matrices derived from Hie theory was tabulated for 
a range of values of the size parameter and refractive Index typical 
of atmospheric particles, integration over the size parameter yielded 
the scattering matrix elanents for a variety of hypothesized particu- 
late size distributions. A least squares curve fitting technique ms 
used to find a best fit from the library data for the experimental 
measurements. This was used as a first guess for a nonlinear Iterative 
Inversion of the size distributions. A real index of 1.50 and an 
imaginary index of -0.005 are representative of the anoothed inversion 
results for the near ground level atmospheric aerosol in Tucs<^. 


I. Introdyctloii 

/ytgulir light seattAring data from tha atmosphare naar ground tevai 
Mara obtained with a polar naphalomatar designed and emstructed at tha 
tbitversity of Ari«>na.^ Tha angular scattering nwasurements made with 
this instruAMnt were processed to acquire the four scattering matrix 
elements due to particulates from a set of four intensity nwasurements 
at each of a nuadter of scattering angle ;. The best data available 
from the instrument were used for analysis of particulate characteris- 
tics. 

*tn order to analyze the experimental results. Hie theory, which 
assumes tosogeneous sfdieres, was applied to develop a theoretical 
particulate scattering representation for comparison. Hie theory 
accurately describes the scatter from the small, typically irregular 
particles found in the atmosphere, however larger particles (in the 
Nie regime) are not particularly well represented. Papers have 

appeared recently on developing a theory to correct for the discrepan- 

3 k 

cies, however further development is necessary before these methods 
can be applied to the extent needed for this study. Therefore, due to 
the difficulties involved in non spherical analysis and the method of 
data acquisition,^ Hie theory was necessarily used for the calculations 
in this paper. 

An inversion of the light scattering data was necessary to obtain 
the best possible representation for the particulates. The application 
of inversion methods for remote sensing has been a region of expanding 
activity since computer development made the necessary manipulations 
tractable. Early techniques included the Phi Ilips-Twomey linear 
inversion^*^ that was applied by Twomey^ to extract vertical temperature 
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prof ties In the sl»M>if^re. This Itn^r method hes been used extensively 

with success In tlM» wnelysts of atmospheric particulates from mu1tlv»ve- 
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length extinction date. 

The Inversifxn of atmospheric aerosol angular scattering data to 

obtain particulate InforiMtltM has typically met with <xi1y marginal 

success. . Vestwater and Cohen^^ felt that the Backus-Gilbert Inversion 

could retrieve size distributions with angular scattering data from 

their theoretical study with mu It {wavelength scattering, fost^^ applied 

this method to imiltlpte angle scattering measurements from water droplets 

but had poor results at sizes below 10pm In radius. Both Post and 

Vestwater and Cohen used narrow size distributions and still had a 

deterioration of results at small sizes* Some success has been achieved 
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in invertii^ bistatic lidar data of atiMspherlc particulates with the 

n 

linear method, but only a very limited data set was available. 

The reasons for using one inversion schmne over another are almost 

as varied as the investigators, however the nonlinear algoritim technique 
lb 

of Twomey has shown promise in retrieving particulate size distribu- 
tions and was chosen on this basis for appllcatltxn to angular scattering 
measurements in this paper. 

II. Data Evaluation Method 

An estimate of the size distribution and refractive index of an 
atmospheric sample is made by comparing the four scattering matrix 
elements measured at various angles with the nephelometer with elements 
produced by theoretical size distributions for various Indices of 
refraction. A library data set on magnetic tape was created using a 
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si^routtn* by Dtv« for scattering by e sphere. Matrix elements were 

recorded for 500 size parameters 10.3(0.2)100.1] at every Integral angle 

[0*(t*)l80*l of scatter for alt eomblnatims of a set of real Indices of 

i^fraetl<m (1*^0, 1.4St 1*50, l.5^» 1.60] and Imaginary Indices of refrae 

tl«i CO.O, -0.003» *0.005» *0.01, *0.03]. (Size parameter, a, is 2n 

times the particle radius divided by the wavelength of Incident light, 
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0.51^5 vn fn this case.) The range of real Indices was chosen to encom* 
pass a regicm from near that of water up past that for silicates. The 
lmaglf»ry Index values vary from ho absorption as for water to a value 
of -0.03 which is near what King* has observed'. 

Subs^uently, these date are Integrated over the size parmneters 
for dunge and two-slope size distributions (Figs. 1*3). The Junge 
size distributions are calculated by setting 
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where N Is the particle number concentration, r is the particle radius, 

C Is a normalization value set to give 100 yg/m^ mass loading, and v 
varies over a typical range frcm 2.0 to 6.0 In 0.2 steps. The two-slope 


'size distributions are calculated by setting 
dH 
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where all combinations of r^ ■ 0.06 jim, r^ - 0.6 and 1.0 wm, ■ 2.0, 
3.0, and 6.0, and V 2 “ 0.0, 1.0, and 1.5 are used. These parameters 
were ctosen to give turnover values between 0.01 and 0.1 yro. These are 
above the values observed by Twomey,^^ but a higher turnover point is 
necessary If any effect were to be observed on the scattering data. 
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GoaparltMi are MKia betw^n these size distributions and real 
data ^ alloMtnf the mss load Inf to vary to five a best least squares 
fit. The quality of the fit i# determined ^ tlw size of H, given by 

M - E tb, - 5d,) (3) 

t * ' 

2 

idiere C is maluated at Ed.b./Ed, for a miniimm H and functicvts as a 

t • ‘ I ‘ 

mss loading .adjusmmt to obtain the best fit. b| is the observed aerosol 
scattering mtrlx elment, and d| Is the correspondli^ theoretical 
mtrix element normalized for 1(KI yg/ro^. Biasing of the data according 
to scattering volume Is also used. Outside the range of the size 
parameters on ta|M, the nuadter concentrations are Inadequate (for any 
realistic size distribution and visible wavelengths of light) to affect 
the observed scatter and are neglected. The tabulated comparisons are 
evaluated to find the best least squares fit with reasonable mss loading 
and to observe any tendencies such as sensitivity to the parameters that 
are varied. 

Due to the similarity of many of the kernels, little Informtlon 
Is gained by using a complete range of angles to obtain a size dlstrlbu' 
tton. The additional time involved In making excessive measurements can 
also be detrimental due to possible changes In the sampled aerosol. 
Therefore, consideration should be given to which angles are most 
critical. Angles vdiere the scattered light Is minimal have more error 
and should be avoided. Also, angles where the scattered radiance changes 
very quickly ore affected more by positioning error in the detector. 

By considering angular scattering measurements mde on monodlsperse 
particles,^ one finds that for larger nonspherical laboratory aerosols. 
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HI* thMfy SBons to hold best for angles less then about 40 degrees. 

^Iler aerosols (as they approach the Rayleleh regime) tend to- follow 
Hie thMry quite well. This leads one to inspect where the larger 
aerosols contribute to the scatter. By looking at Figs. 4 and one 
observes that the difference due to the targe aerosols is limited 
iMiinly to the forward few degrees. This not only inq>1tes that Hie theory 
should hold better for a typical aerosol size distribution than for 
tingle large aerosol studies, but that if one desires information conr 
tent from the larger aerosols, measurements must be made in the forward 
direction or little Information above I ym Is obtained for typical size 
distributions. 

III. inversion Technique 

Continuing with the next step, the Inversion method is considered. 
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A nonlinear algorithm was developed by Chahine which essentially 

4 

assumed delta functions for kernels but acquired inherent instabilities 

due to increased high frequency content when measurements were numerous. 

14 

Chahine's algorithm was modified by Twomey to include the entire 
nonzero region of the kernel. This eliminated the detrimental factor of 
superfluous data and. In fact, caused the inversion to improve with 
additional data due to an effective decrease in measurement error. 

The nonlinear Inversion has also shown an ability to cope with measure- 
ment errors, which greatly strengthens Its position in application to 
the aerosol size distribution problem. The iterative algorithm Is 
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where f (r) ts the Inttiel guess size distribution, k(r,s) Is the kernel 
value (the theoretical scattering matrix element for a single particle), 
g(s) ts the actual measured scattering matrix element for the collection 
of particles, and f'(r) Is the modified size distribution. Variable r 
refers to the particle radius or size parameter, and variable s refers 
to a specific matrix element measurement. Examination of the kernels 
shows that a lot of fine structure typically occurs (Figs. 6-9), particu- 
larly near backscatter. Whereas this might be expected to be beneficial 
for fine resolution, in practice this structure is too high a frequency to 
be effective in Improving the Inversion accuracy. Since atmospheric aero- 

j 

sol size distributions do not seon to have these wild oscillations and 
neither do the observed scattering measurements, the fine structure would 
not seem necessary to resolve that data even if it were effectively usable. 
In fact. It might be desirable to use smoothing of the kernel to assist In 
obtaining a stable solution. The power, spectra of the kernels also stow 
that the middle frequencies are often deficient (Fig. 10), and occasionally 
even low frequencies are absent (Fig. II). This Is a strong negative factor 
In the application of scattering kernels to Inversion techniques. 

Simple quadrature Is used for the Integral with the kernels being 
read from magnetic tape. Each data value is successively Iterated once 
through all the particle sizes on tape, modifying the size distribution 
according to the kernel's weighting effect. The weighting is scaled to 
less than or equal to one by dividing by the maximum kernel value for a 
particular angle and matrix element. After each unknown in the set has 
been determined from the first Iteration, the process is repeated until 
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the final distribution Is obtained. Although the inversion lends Itself 
easily to programming, care is still required in its application. 

IV. Inversion of Theoretical Data 

Initial runs of the Inversion program were made on data generated 
from Mie theory to establish the accuracy of the Inversion with scattering 
kernels. The nephelometer measures the radiance of tight scatter which 
is a function of the particles' scattering cross sections times their 
concentrations. It was necessary to weight the scattering kernels 
according to an initial, first guess size distribution to obtain reasonable 
results. Otherwise, there was a strong tendency for the inversion to 
adjust the large particle concentrations to the point of instability. 

The runs were made using only the M 2 and Mj elements^ from five forward 
angles and two backward angles. These were chosen to maximize information 
content with a minimum of data. A special problem occurs in applying 
Inversions to the $ 2 | and 02 ^ elements as it is possible for the theoreti- 
cal and measured values to be of opposite sign due to errors In the 
measurements or in the first guess size distribution. This would imply 
a negative particle concentration that is not allowed. Runs were made 
with theoretical data from Junge distributions. A v of two, index of 
1.54- .0051, and mass loading of 38 pg/m^ were used for a first guess, 
as these values produced a close fit for one of the real aerosol runs. 

A method of overrelaxation was settled upon as the best technique for 
applying the algorithm, it has the format 

f(r) - {I + m’'} f°(r) (5) 
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where 


N.Ur—aisl I 

J2f (r)k(r,s)dr „ 

The absolute value of M Is raised to the power R, and that quantity takes 
the same positive or negative sign as M. 

Not only did overrelaxation with values of R less than one speed 
up convergence, but it Improved the results greatly (Fig. 12). Too large 
of an overrelaxation, however, caused oscillations. A value of 0.7 for R 
produced the best stability and convergence although 0.5 gave the fastest 
convergence. Excessive Iterations are not only costly, but tend to pro- 
duce a more highly structured, atypical size distribution. This Is 
avoided by terminating the Iterative process after successive Iterations 
with less than 0.2 percent Improvement In error. 

The Inversion program was run with theoretical data to observe the 
effect of various size distribution first guesses on the results of the 
Inversion. Large differences between the actual and Initial guess mass 
loading were difficult for the Inversion to handle If no Internal mass 
loading adjustment Is Included in the program (Fig. 13). This serves 
to point out graphically the size range of Information content of the 
data and the kind of structure that can occur due to the oscillatory 
nature of the kernels. The response In the region of Information 
content Is adequate to Indicate the proper correction necessary (I.e. 
higher or lower) In the mass loading. Differences between the actual 
and initial guess v values (fig. |l>) have much less effect on the 
Inversion. Examination of the region of convergence does show that the 
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scattering is mainly sensitive to particles in the 0.2 to 2 ym range 
and subsequently, this is the region vdiera results are applicable. 

Random error vies added to the theoretical data to ^serve its 
effect on the inversion. Neither l| nor II percent error had a signifi- 
cant effect on the inverted site distrilMttion. This result Is essential 
to obtaining realistic inverti<ms with experimental data. Tlw convergence 
limit on the error (Fig. 1$) was, in fact, indicative of the percentage 
error in the data; however, more study of this point is necessary for 
verification. 

The linear inversion method was also applied to this problem initially. 
However, it could not invert the data unless the error level was I percent 
or less. This is an unrealistically low value, especiatly since Mie theory 
alone can account for more than 1 percent 'error. Therefore, the linear 
method was dropped. 

V. Experimental Results 

After checking the ability of the inversion to reproduce theoretical 
data, experimental data were analyzed. 

A. Curve Fitting 

By comparing the experimental data with the theoretical library 
data, a best fit was obtained. Typically, about 20 scattering angles 
with four matrix elements at each angle were used. The fit was weighted 
by the cosecant of the scattering angle to allow more bias for larger 
scattering volumes. After checking the data fit with matrix elements 
produced by both Junge and two-slope size distributions and using trun- 
cated data sets that excluded measurements of smaller magnitude, an 
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MtlMts for ^ MroMl s!m distribution was obtain^. The nepheionwtar 
runt for iiarch $th Mere averaged and gave a Junge best fit with 
m " 1.50 • .0031* V " 2.1* and iMtt loading ■ 38 vg/M^. The ^-siope 

best fit for the taise data was m ■ 1.50 - .0031* vj • 2.0* V 2 ■ 0* 

r^ ■ .04 |M}» r^ ■ 1.0 )m* and mast loading ■ 80 ug/n • 

The mass loading is not extrmnely critical as the largest particles 

dominate this value while they have much less effect on the actual 

light scatter on tdilch the i^surmsents are bas^. Another set of 

runs on Harch 10th was average to yield a Jtmge best fit with 

m * 1.47 * .0051* V ■ 2.0, and mass loading " 47 vg/m and two-slope fit 

with m ■ 1.50 - .0041, vj • 2.0, V 2 • 0,0, r^ ■ .04 »«», r^ 1.0 >«d, 

3 

and moss loading ■ 100 pg/m . Figs. 16 end 17 show typical graphs of 
the experimentally measured matrix eleiMnts plotted in comparison with 
the theoretical data produced by the corresponding best fit size dlstrlbu- 
ti(W). As expected, tiM curves match closely near the forward direction 
which is where Hte theory is believed to hold best and where the strongest 
weighting is placed on the least squares fit. The two-slope end Junge 
distributions vdiich gave the test fits are of similar form over the size 
range of information content. Therefore, the simpler dunge distribution 
was chosen for coi^lning the March data which gave an overall aerosol 

3 

characterizatifw) of m - 1.49 - .004t, v > 2.0, and mass loading • 40 pg/m . 

The strongest sensitivity for the ranges of parmneters under con- 
sideration was observed to be the size distribution slo^^.., next In 
i^x>rtance tMS the imaginnry Index, and the least sensitive was the 
real Index. 
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I. Inversion of Site Distributions 

Further Inprovenent In the size distribution estimete is attmpted 

by using the llbrery ^ta closest to the curve fitting results as a 

* 

first guess for the Inversion. Fairly consistrnt results are obtained 

by the Inversion (as shown by Ftg. 18 for the March 9th data) even if the 

Initial V value is varl^ above or below the curve fitting results. 

The convergence of the Inversiwi Is shoturn by the error smoothly 

approaching minimum values as the number of iterations Increases (Ftg. 19) • 

Even when the initial guess is greatly in error from the data, the RMS 

error iterates down to the 15 percent range which is representative for 

all the runs. Attaints were made to improve the minimum RMS error of 

the Iterated Inverstwi by using different indices of refraction; however, 

this exercise just verified the choices of the curve matching technique. 

Frcd>able causes of this large of a convergence limit (If It is truly 

Indicative of the experimental error) are covered In a preceding paper.* 

OtiMr Inverted size distributions with various initial guesses are 

shown in Figs. 20 and 21. Twelve inversions were averaged and smoothed 

to obtain a representative Inverted size distribution (Fig. 22). The 

results are most closely modeled by a Junge size distribution with 

s 

m ■ 1.50 - .0351, V ■ 1.8, end mass loading of about 80 yg/m . Ti^ curve 
is purposely truncated so that only the regicwi of sensitivity is shown. 

A maximjm aerosol number concentration (or turnover point) of the size 
distribution is not observed since the sensitivity of the kernel drops 
off sharply below 0.2 ym tdille typical turnover points occur near 0.01 ym 
for the ground level Tucson aerosol. 
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VI. Conclusions and Further Study 

The technique developed in this paper has yielded estimates of 
atmospheric aerosol characteristics — vis.» size distributions and real 
and imaginary indices of refraction— from measurements of matrix 
scattering elements at various angles. The results are reasonable in 
comparison with other work in this field, and the size distributions 
match quite well near 0.1 ym with nuclepore filter measurements^^ that 
were made at the same location but are sensitive to particles from 
0.1 ym <Mi dcxfn. To achieve a characterization of the ground level 
Tucson aerosol, measurements should be made routinely over an extended 
period. 

Improvement in the stability of the inversion technique might be 
achieved by smoothing the kernels to remove the higher frequency informa- 
tion. Further sophistication could be achieved by expanding the inversion 
routine to include fitting the inversion results to a smooth analytic 
function and using this as a new first guess. 

A detailed study of the information content would be of special 
interest. Initial work In this area has shown that the scattering data 
used in this study are in the region of optimum information content. 

This could lead to an instrunent with a minimum number of fixed detectors 
set at carefully chosen angles which would eliminate the need to move 
the detector and speed the measurement time. 

The author would like to thank Dr. Benjamin M. Herman and Dr. Walter H. 
Evans for their assistance on this project. This research was funded by 
the Office of Naval Research under Grant #NOO0l4-75-C-O2O8, and computer 
time was furnished by the National Center for Atmospheric Research which 
is sponsored by the National Science Foundation, under project #35021004. 
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Figure Captions 


Fig. 1. Junge size distribution curves. 

Fig. 2. Two-slope size distribution curves. 

Fig. 3. Scattering element H2 integrated over Junge size distributions. 
Fig. k. Scattering element S2^ integrated over a Junge size distribution 
for two particle size ranges. 

Fig. 5. Scattering element M2 integrated over a Junge size distribution 
for two particle size ranges. 

Fig. 6. Weighted scattering element M2 for single particles. 

Fig. 7. Weighted scattering element M^ for single particles. 

Fig. 8. Scattering element M^ for single particles. 

Fig. 9. Scattering element M2 for single particles. 

Fig. 10. Power spectrum of M2* 

Fig. 11. Power spectrum of D2J. 

Fig. 12. Convergence of iterative inversion for theoretical data with 
no error. 

Fig. 13* Theoretical size distribution inversions for various mass 
loading initial guesses. 

Fig. 14. Theoretical size distribution inversions for various Junge 
slope initial guesses. 

Fig. 15. Convergence of iterative inversion for theoretical data with 
1 1 percent error. 

Fig. 16. M2 scattering matrix element. 

Fig. 17. $21 scattering matrix element. 

Fig. IS. Size distribution inversions for March 9th experimental data. 


Ftg. 19* Convergence of Iterative Inversion for experimental data with 
various mass loading Initial guesses. 

Fig. 20. Size distribution Inversions for March 4th experimental data. 
Fig. 21. Size distribution Inversions for March lOth experimental data. 
Fig. 22. Average of inverted size distributions. 
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Fig. 4 . Scattering element S21 integrated over a 
Junge size distribution for two particle sizi 
ranges. 
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Fig. 6. Weighted scattering element M 2 
for single particles. 
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Fig. 13. Theoretical size distribution inversions 
for various mass loading initial guesses. 
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Fig. 14. Theoretical size distribution inversions for 
various Junge slope initial guesses. 
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Fig. 16. M 2 scattering matrix element. 
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Fig. 17. $21 scattering matrix element . 
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Fig. 18. Size distribution inversions for March 91 
experimental data. 
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Fig. 20. Size distribution inversions for March 4tl 
experimental data. 
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Fig. 21. Size distribution inversions for March 1C 
experimental data. 
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